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ABSTRACT 

We present a new public archive of light-motion curves in Sloan Digital Sky Survey (SDSS) 
Stripe 82, covering 99° in right ascension from a — 20.7^ to 3.3*^ and spanning 2?52 in 
declination from S = — 1?26 to 1?26, for a total sky area of ^249 deg^. Stripe 82 has been 
repeatedly monitored in the u, g, r, i and z bands over a seven-year baseUne. Objects are 
cross-matched between runs, taking into account the effects of any proper motion. The re- 
sulting catalogue contains almost 4 million light-motion curves of stellar objects and galax- 
ies. The photometry are recalibrated to correct for varying photometric zeropoints, achieving 
~20 mmag and '^30 mmag root-mean-square (RMS) accuracy down to 18 mag in the g, r, i 
and z bands for point sources and extended sources, respectively. The astrometry are recali- 
brated to correct for inherent systematic errors in the SDSS astrometric solutions, achieving 
^^32 mas and '^35 mas RMS accuracy down to 18 mag for point sources and extended sources, 
respectively. 

For each light-motion curve, 229 photometric and astrometric quantities are derived and 
stored in a higher-level catalogue. On the photometric side, these include mean exponential 
and PSF magnitudes along with uncertainties, RMS scatter, per degree of freedom, various 
magnitude distribution percentiles, object type (stellar or galaxy), and eclipse. Stetson and 
Vidrih variability indices. On the astrometric side, these quantities include mean positions, 
proper motions as well as their uncertainties and per degree of freedom. The here presented 
light-motion curve catalogue is complete down to ?■ ^ 21.5 and is at present the deepest 
large-area photometric and astrometric variability catalogue available. 

Key words: catalogues - stars: photometry, astrometry, variables - Galaxy: stellar content - 
galaxies: photometry 



1 INTRODUCTION 

One meaning of the verb "to vary" is to change in amount or level, 
especially from one occasion to another. In astronomy, there are 
many types of temporal variability. Stars may change in brightness, 
in which case they are termed "variables", or they may change po- 
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sition on the sky, in which case they have a proper motion. Over a 
sufficiently long period of time, the shape of constellations change. 
Galaxies exhibit variability across the electromagnetic spectrum 
since their emission is made up of the radiation from billions of 
sources, although their most obvious source of variation comes 
from active galactic nuclei or supemovae. Transient events, such as 
gamma-ray bursts, or microlensing events, are instrinsically vari- 
able. Solar system objects from planets to asteroids drift slowly 
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Table 1. The list of SDSS imaging runs included in the LMCC. All runs were processed with version 40 of the SDSS frames pipeline except those runs 
marked with an asterisk, which were processed with version 41. Runs 4203 and 5823 (underlined) are the reference runs used for the astrometric calibrations 
(see Section 2.4). 



Year Month North Strip Runs South Strip Runs 



1998 


Sep 


94 








125 






1999 


Oct 


1033 








1056 






2000 


Sep 


1752 














2000 


Oct 










1755 






2000 


Nov 










1894 






2001 


Jun 


2385 














2001 


Sep 


2570 


2578 


2589 




2579 


2583 


2585 


2001 


Oct 


2649 


2650 


2659 


2662, 2677 








2001 


Nov 


2700 


2708 


2728 


2738 


2709 






2001 


Dec 


2768 


2820 












2002 


Jan 


2855 


2861 


2873 




2886 






2002 


Sep 










3325' 




2002 


Oct 


3362 


3384 


3437 




3355 


3360* , 3388, 3427, 3430. 3434. 3438 


2002 


Nov 


3461 








3460 


3465 




2003 


Sep 


4128 


4153 


4157 




4136 


4145 




2003 


Oct 


4184 


4188 


4198 


4207 


4187 


4192 


4203 


2003 


Nov 


4253 








4247 


4263 


4288 


2004 


Aug 


4797 














2004 


Sep 


4858 














2004 


Oct 


4868 


4874 


4895 


4905,4917 








2004 


Nov 


4933 


4948 






4930 






2004 


Dec 










5042 


5052 




2005 


Sep 


5566 


5603 


5610 


5622, 5633, 5642, 5658 


5582 


5597 


5607, 5619, 5628, 5637, 5646, 5666 


2005 


Oct 


5709 


5719 


5731 


5744, 5759, 5765, 5770, 5777, 5781, 5792, 5800 


5675 


5681 


5713, 5729, 5730, 5732, 5745, 5760, 5763, 5771, 5776, 5782, 5786, 5797 


2005 


Nov 


5813 


5823 


5842 


5865, 5866, 5872, 5878, 5898, 5902, 5918 


5807 


5820 


5836, 5847, 5853, 5870, 5871, 5882, 5889, 5895, 5905, 5924 



across the sky, waxing and waning on various timescales. In fact, 
at some level, everything in the sky is variable. 

The introduction of CCD detectors to astronomy greatly en- 
hanced the ability to conduct variability surveys. The extension of 
CCD cameras to mosaic and wide-field formats along with the ex- 
ponential progression of computing power have allowed the subse- 
quent development of more ambitious surveys reaching to deeper 
magnitudes, higher cadences and larger sky area s. For more de- 
tails we direct the reader to iBecker et al.l j2004h who present a 
clear summary of modem variability surveys. In this work we con- 
centrate on optical photometric and astrometric variability (hence 
"light and motion") over a ~249 deg^ patch of sky. 

Large sky sur veys such as the Sloan Digital Sky Sur- 
vey (SDSS; IVork et al. 2000) have in many ways revolutionised 
our knowledge of the Universe. SDSS has imaged approximately 
a quarter of the sky in five photometric wave bands. The exploita- 
tion of this impressive dataset has resulted in hundreds of pub- 
lications covering a wide range of astronomical topics, from the 
structure of the Milky Way to the mapping of a large fraction of 
the Universe. The bulk of this data, however, contain only single 
measurements of objects from the north Galactic cap with no infor- 
mation on possible photometric varia bility or astrometri c motion. 
Substantial efforts have be en made bv lMunn et al.l ( l2004h (see also 
iGould & KoUmeieJ ( l2004f) ) to measure proper motions by match- 
ing S DSS data from the n orth Galactic cap with the USNO-B cata- 
logue jMonet et al.l2003l) . The resultant proper motion catalogue is 
90% complete down to <; = 19.7 with the magnitude limit being 
set by the USNO-B catalogue faint magnitude limits. 

One of the primary goals of the SDS S is the study of the vari- 
able sky jAdelman-McCarthv et al.ll2007h of which our knowledge 
is still very incomplete ( |Paczviiskil2OO0l) . To this end, the SDSS has 



repeatedly imaged a 300 square degree area, the so called Stripe 82, 
during the later half of each year since 19 98. In 2005, the SDSS- 
II Supernova Survey dFrieman et al.ll2008h started with the aim of 
detecting Type-I supemovae in Stripe 82, greatly improving the ca- 
dence of measurements within the stripe. By averaging a subset 
of the repeated ob servations of unresolved sources in Stripe 82, 
llvezic et al. I ( l2007h built a standard star catalogue containing ~1 
million nonvariable sources with r band magnitudes in the range 
14-22, by far the deepest and most numerous set of photomet- 
ric stan dards available. U sing these same multi-epoch photomet- 
ric data. lSesar at al.l ( l2007h analysed the photometric variability for 
~1.4 million unresolved sources in the stripe, drawing interesting 
conclusions on the spatial distribution of RR Lyrae stars and the 
variability of quasars. 

Here we present a new public archive of light-motion curves 
in SDSS Stripe 82. The archive has been constructed from the set of 
high-precision multi-epoch photometric and astrometric measure- 
ments made in the stripe since the first SDSS runs in 1998 until 
the end of 2005. In constructing the catalogue, we only use mea- 
surements of objects that are cleanly detected in individual SDSS 
runs. The catalogue contains almost 4 million objects, galaxies and 
stars, and is complete down to magnitude 21.5 in it, g, r and i, and 
to magnitude 20.5 in z. Each object has its proper motion calcu- 
lated based only on the multi-epoch SDSS J2000 astrometric mea- 
surements. The catalogue reaches almost two magnitudes deeper 
than the SDSS/USNO-B catalogue, making it the deepest large- 
area photometric and astrometric catalogue available. 

The catalogue comes in two flavours, the Light-Motion Curve 
Catalogue (LMCC), which contains the set of individual light- 
motion curves, where measured quantities for each object are listed 
as a function of wave band and epoch, and the Higher-Level Cat- 
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alogue (HLC), which presents a set of derived quantities for each 
light-motion curve. For many purposes, it is more convenient to 
work with the HLC, especially for selecting subsets of interesting 
objects. The construction, calibration and format of the LMCC are 
discussed in Section 2, and the HLC is described in Section 3. Be- 
tween the two subcatalogues, there is all the necessary information 
available to explore the photometric and astrometric variability of 
~249 square degrees of equatorial sky. In Section 4, we investi- 
gate the quality of the photometric and astrometric properties of our 
catalogues by comparing them against suitable external catalogues, 
and we analyse the behaviour of our proper motion uncertainties. 



2 THE LIGHT-MOTION CURVE CATALOGUE 
2.1 Stripe 82 Data 

The SDSS photometric camera is mounted on a 2.5m dedicated 
telescope at the Apache Point Observatory, New Mexico. It consists 
of a photometric array of 30 SITe/Tektronix CCDs, each of size 
2048x2048 pixels, arranged in the focal plane of the telescope in 
six columns of five chips each tGurm et alJl998l : lGunn et al .120061) 
with a space of approximately one chip width between columns. 
Each row of six chips is positioned behind a different filter so that 
SDSS imaging data is produced in five wave bands, namely, u, g, r, 
i and z (Fiikugita et al. 1996; S mith at a l. 2002.) . The camera oper- 
ates in time-delay-and-integrate (TDI) drift-scan mode at the side- 
real rate and the chip arrangement is such that two scans cover a 
filled stripe 2? 54 wide, with '^l' overlap between chip columns in 
the two scans. In addition, the camera contains an array of 24 CCDs 
with 400x2048 pixels which enable observations of bright astro- 
metric reference stars for subsequent astrometry and focus moni- 
toring. 

The images are automatically processed t hrough spe- 
cialised pipelines i Lupton. G unn & Szalav 1999; Lupton et al. 



1200 ll ; iHogg et al.ll200ll ; IStoughton et al. .2002 ; Jvezic et aL,2004.) 
producing corrected images, object catalogues, astrometric solu- 
tions, calibrated fluxes and many other data products. The object 
catalogues, which include the calibrated photo metry an d a strom- 
etry, are stored in FITS binary table format JWells et al]|l98ll ; 
ICotton et all l 19951 ; iHaiiisch et al.ll200lh and referred to as "tsObj" 
files. It is these object catalogues that we have used to construct the 
LMCC. 

The SDSS Stripe 82 is defined as the region spanning 8 hours 
in right ascension (RA) from a = 20"^ to 4^ and 2? 5 in declination 
(Dec) from 5 — — 1?25 to 1?25. The stripe consists of two scan 
regions referred to as the north and south strips. Both the north and 
south strips have been repeatedly imaged from 1998 to 2005, be- 
tween June and December of each year, with 62 of the 134 imaging 
runs obtained in 2005 alone (this large sampling rate was produced 
by the start of the SDSS-II Supernova Survey). A specific imag- 
ing run may cover all of one strip or some fraction of the area, and 
images of the same patch of sky are never taken more than once 
per night. Hence the exact temporal coverage and cadence of the 
light-motion curves in the catalogue are strong functions of celes- 
tial position. In Table[T] we list the SDSS imaging runs included in 
the LMCC organised by the month and year of observation. Not all 
scans of Stripe 82 were inclu ded in our analysis due to failures in 
the SDSS frames pipeline ( Lupton et alj|200lh when processing 
a run or failure of our calibration routines to produce photometric 
zeropoints (see Section 2.2). 



0.04 
0.02 

I -0.02 
o 

^ -0.04 

CD 
in 

5 -0.06 
^ -0.08 

CD 
CD 

-0.12 
-0.14 
-0.16, 



150 200 250 300 350 400 450 

Field Number 

Figure 1. Plot of the median fractional flux offset of the reference stars 
relative to the reference catalogue as a function of field number for runs 94 
(solid line) and 5853 (dashed line). For both runs we show the offsets for 
camera column 1 and the r band. Run 94 is a typical pre-2005 run observed 
under photometric conditions and run 5853 is a typical 2005 supernova run 
observed under non-photometric conditions. 



2.2 Further Photometric Calibrations 

The Stripe 82 data set includes 62 "standard" SDSS imaging runs 
which were observed under photometric conditions and which 
were photometrically calibrated using the standard SDSS pipelines 
(Tucker et al. 2006). We use these standard runs to construct a ref- 
erence catalogue of bright star fluxes, from which we can both im- 
prove the photometric calibrations of the standard runs, as well as 
derive photometric calibrations for the Stripe 82 supernova imag- 
ing runs from 2005, which were generally observed under non- 
photometric conditions. 

To construct the reference catalogue, we start with a set of 
bright, unsaturated stars, with 14 < r < 18, taken from a set 
of high quality runs acquired over an interval of less than twelve 
months (2659, 2662, 2738, 2583, 3325, 3388). We then match the 
individual detections of these stars in each of the 62 standard runs, 
using a matching radius of 1 arcsec. On average, there are 10 inde- 
pendent measurements of each star among the standard runs, and 
we only include in the reference catalogue stars with 5 or more 
measurements. We then compute the unweighted mean of the inde- 
pendent flux measurements of each star and adopt that mean flux 
in our reference catalogue. Note that we specifically use the fluxes 
measured in the so-called SDSS "aperture 7", which has a radius 
of 7.43 arcsec; this aperture is typically adopted in the SDSS as a 
reference aperture appropriate for isolated bright star photometry. 

These reference catalogue stars are used to re-calibrate the 
standard runs, as well as to calibrate the 2005 supernova data. For 
the supernova runs, we first adopt a sensible but arbitrary zeropoint 
for the purposes of generating the initial tsObj files using standard 
SDSS pipeline tools. We then match the object detections in each 
run to the reference catalogue, and compute the median fractional 
flux offset of the reference stars in the individual run relative to the 
reference catalogu^ These median offsets are computed in two it- 



^ The u band images possess a significantly poorer signal-to-noise ratio 
than the other bands and the use of reference catalogue stars that have a 
magnitude fainter than 18 before calibration (due to higher than usual ex- 
tinction) degrades the determination of the derived flux offsets. The flux 
offsets for the u band have therefore been determined using only the ref- 
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Table 2. The list of extra constraints that all need to be satisfied in at least one wave band in order for an object record to be included in the LMCC. 



Tag Name In Relation Flag Value Flag Name 

tsObj File (Hexadecimal Bit) 



OBJCJLAGS 

FLAGS 

FLAGS 


AND 
AND 
AND 


0x4 
0x10000000 
0x20 


FALSE 
TRUE 
FALSE 


FLAGS 


AND 


0x80 


FALSE 


FLAGS 
FLAGS 


AND 
AND 


0x40000 
0x80000 


FALSE 
FALSE 


FLAGS 


AND 


0x400000 


FALSE 


FLAGS2 


AND 


0x100 


FALSE 


FLAGS2 
FLAGS2 
FLAGS2 
FLAGS2 


AND 
AND 
AND 
AND 


0x800 
0x1000 
0x4000 
0x8000 


FALSE 
FALSE 
FALSE 
FALSE 



erations. We first calculate the median fractional offset for each run 
in bins of 0.0208° in Dec, i.e., 120 bins over the width of Stripe 82, 
or about 10 bins per CCD width. This exercise is designed to correct 
flatfielding errors for a given run. Note that these errors would only 
depend on Dec because the SDSS employs a drift-scan camera, and 
the scan direction for Stripe 82 is in the RA direction. After cor- 
recting for the declination-dependent offsets, we then re-compute 
the median fractional flux offsets for each field along a given run 
(each SDSS field is 0.15° long in RA). This additional field-by- 
field offset corrects for any temporal variations in the photometric 
zeropoint of a given run, which are due to transparency/extinction 
changes over the course of a nominally photometric night. 

Flux offsets for a certain wave band were only applied to ob- 
jects assigned to a bin with at least 9 reference catalogue stars in 
order to guarantee the accuracy of the derived flux offset. In prac- 
tice, this extra restriction only affects the photometry of objects in 
the u band, and in other wave bands when the atmospheric trans- 
parency is low. We use a photometric calibration tag (see Table[3} 
to monitor whether or not a flux offset has been applied to the pho- 
tometry of a particular object at a certain epoch in a specific wave 
band. The final tsObj files used for subsequent analyses therefore 
have both these declination-dependent and field-dependent flux off- 
sets removed for most object records. We find that for the standard 
SDSS runs, the fractional flux offset corrections, which we refer to 
as photometric zeropoints, are about 1 — 2%, which sets the typi- 
cal scale of these residual errors in the standard SDSS calibration 
procedures. Figure[T]shows the fractional flux offset corrections as 
a function of SDSS field number (an arbitrary coordinate along RA 
assigned to image sections from the same run) for a typical photo- 
metric run (94) and a typical non-photometric run (5853). 



erence catalogue stars with an uncalibrated magnitude brighter than 18 in 
each run. 



Description 



Reject objects too close to the edge of the image 

Accept only objects detected in the unbinned image 

Reject objects where the given centre is the position of the peak 

pixel, rather than based on the maximum likelihood estimator 

Reject objects that are too small or too close to the edge to 

estimate a radial profile 

Reject objects with one or more saturated pixels 

Reject objects with pixels that were not checked to see whether 

they included a local peak, such as the cores of saturated stars 

Reject objects with a sky level so badly determined that the highest 

pixel in an object is very negative, far more so than a mere 

non-detection 

Reject objects containing interpolated pixels that have too few 

good pixels to form a reliable estimate of the flux eiTor 

Reject objects with a centre close to at least one saturated pixel 

Reject objects with a centre close to at least one interpolated pixel 

Reject child objects with no detected peak 

Reject objects with more than 20% of the PSF flux from 

interpolated pixels 



2.3 Catalogue Construction 

The object catalogues (tsObj files) contain quality and type flags 
for each object record to aid in the selection of "good" measure- 
ments and specific data samples. In the LMCC, we only accept 
object records classified as galaxies/non-PSF-like objects (tsObj 
file tag OBJC.TYPE = 3) or stars/PSF-like objects (tsObj file tag 
OB JC.TYPE = 6), a nd the object mus t have no child objects (tsObj 
file tag NCHILD = 0:^ Stoughton et aL Il l2002h ). We then require that 
an object record satisfies all of a set of constraints in at least one 
wave band. The first of these constraints is that a photometric ze- 
ropoint, calculated using the method described in Section 2.2, has 
been applied to the object record, and that the object record has an 
uncalibrated PSF magnitude (tsObj file tag PSFCOUNTS) brighter 
than 21.5 for the bands u, g, r and i, or brighter than 20.5 for the 
2 band. These limits were chosen to ensure that any photometric 
measurement in the LMCC has a signal-to-noise ratio of at least 5 
in at least one wave band. In Table (2] we list the remaining set of 
constraints to be satisfied in at least one wave band in order for an 
object record to be included in the LMCC. 

We apply one final constraint on the quality of an object record 
in order to avoid the inclusion of cosmic ray events in our cata- 
logue. If an object record satisfies all of the above constraints in 
one wave band only, then it is accepted only if the tsObj file tag 
FLAGS 2 for that wave band does not contain the hexadecimal bit 
0x1000000 (flag name MAYBE.CR), the presence of which indi- 
cates that the object is possibly a cosmic ray. 

In order to construct the light-motion curves, we processed 
each run in turn, starting with the 2005 runs which were closely 
spaced in time. For each object record in the current run satisfying 
our quality and type constraints, we used the following algorithm 
to process the record: 

(i) We define a subset of objects with light-motion curves from 
the current catalogue that have mean positions inside a 1' box cen- 
tred on the position of the current object record. 

(ii) We calculate an expected position at the epoch of the cur- 
rent object record for each object in the subset using one of two 
methods depending on the number of epochs in the corresponding 
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Figure 2. Top panel: Greyscale image showing the maximum number of epochs in the light-motion curves as a function of object position. Bottom panel: 
Maximum number of epochs in the light-motion curves as a function of RA for a 0?01 wide slice through the greyscale image centred at 5 = — 0?1 
(continuous line; slice through the North strip) and for a similar slice centred at 5 = 0? 1 (dashed line; slice through the South strip). 



1000 F 



100 



00 



1 



1 


1 


1 , , , 1 , , , 1 _ 


; 10000 


Stars 


■ ■ . / 

. .■■ ..." . .■ 








1 


1 





14 


1 6 


18 20 22 




Mean 


r PSF Magnitude 



1000 r 



10000 Galaxies 



100 r 




1 r 



14 16 18 20 

Mean r Exponential Magnitude 



22 



Figure 3. Left: Plot of RMS PSF magnitude deviation versus mean PSF magnitude for 10000 random PSF-like objects (stars) with at least 20 good r magnitude 
measurements. Right: Plot of RMS exponential magnitude deviation versus mean exponential magnitude for 10000 random non-PSF-like objects (galaxies) 
with at least 20 good r magnitude measurements. Both panels: Plotted functions (continuous lines) are all of the form /(m) = A + B exp (C (m — 18)) 
where m denotes magnitude and A, B and C are fitted parameters. For the left hand panel. A, B and C have values 17.3 mmag, 2.48 mmag and 1.00, 
respectively. For the right hand panel. A, B and C have values 30.8 mmag, 5.67 mmag and 1.01, respectively. 



light-motion curve. If an object has a light-motion cui^e with six 
epochs or less, then the mean position is used for the expected po- 
sition. Otherwise a mean position and proper motion are fitted to 
the light-motion curve and used to calculate the expected position 
of the object at the epoch of the current object record. 

(iii) From the expected positions of the subset of objects, we 



find the closest object to the position of the current object recorcQ 
If the closest object lies within 0.7", then the current object record 



^ Note that an object record contains a single datum for the celestial coor- 
dinates, calculated from the astrometric solution for the SDSS CCD camera 
at the current epoch. 
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is appended to the light-motion curve of the closest object, other- 
wise a new light-motion curve is created containing only the current 
object record. 

Since each run contains at most one measurement of any object, the 
above algorithm can be performed in parallel for all objects from 
one run. 

We choose to include both the PSF magnitude and 
exponential-profile magnitude jStoughton et a l. 2002) in a light- 
motion curve as a measure of an object's brightness in each wave 
band at each epoch. The PSF magnitude is the optimal measure of 
the brightness of a point-source object, and hence it is suitable for 
studying stars and quasars. Photometry of extended objects, such 
as galaxies, may be performed in a variety of ways, including fit- 
ting an exponential profile to the object image. The advantage of 
including the exponential magnitude as opposed to any of the other 
available profile magnitudes is that the difference between the PSF 
and exponential magnitudes, referred to as a concentration in dex, 
may b e used as a continuous object-type classifier ( Scrant on et al.l 
l2002h . independent of the more restrictive binary SDSS classifica- 
tion. 

After processing all Stripe 82 data, the LMCC was trimmed to 
only those objects that have a mean position in the range a = 20.7'' 
to 3.3*" and 5 = -1?26 to 1?26. This was desirable because the 
temporal coverage is too sparse outside these limits. The total area 
of sky covered by the LMCC is therefore ~249 deg^. The LMCC 
was also searched for photometric outliers using a 3(7-clip algo- 
rithm for outlier identification, and we found clear groups of out- 
liers clustered at specific HJDs. The tight clustering in time indi- 
cates that these groups of outliers are simply due to bad epochs (and 
not due to some astrophysical process, like an eclipse or flare), and 
hence we removed the corresponding data points, which amounted 
to ^^1% of all epochs. 

In the top panel of Figure |2l we plot a grey scale image show- 
ing the maximum number of epochs in the light-motion curves as 
a function of object position. One may clearly see that the light- 
motion curves for objects in the overlap regions between the north 
and south strips contain approximately twice as many epochs as 
the light-motion curves for objects elsewhere. In the bottom panel, 
we plot a one-dimensional slice through the greyscale image for 
5 = — 0? 1 and 5 — 0? 1 to further illustrate the dependence of 
the number of light-motion curve epochs on RA. 

In Figure [3] we plot sample photometric RMS diagrams for 
the r band. The left hand panel shows the RMS PSF magni- 
tude deviation versus mean PSF magnitude for 10000 random 
PSF-like objects (MEAN.OB JECT.TYPE = 6; see Section 3.1) 
that have light-motion curves with at least 20 good r mag- 
nitude measurements. PSF-like objects are mainly stars with 
some contamination by quasars. Similarly, the right hand panel 
shows the RMS exponential magnitude deviation versus mean 
exponential magnitude for 10000 random non-PSF-like objects 
(MEAN_OBJECT_TYPE = 3) that have light-motion curves with 
at least 20 good r magnitude measurements. Non-PSF-like objects 
are mainly galaxies. 

Overplotted on each RMS diagram is an empirical fit 
to the data using an exponential function of the form 
/(m) = A + B exp (C (m — 18)) where m denotes 
magnitude and A, B and C are fitted parameters (whose values are 
reported in the caption of Figure [Sj. The fitting was done using an 
iterative 3cr-clip algorithm (see Vidrih index under Section 3.1). It 
is clear from these diagrams that for the r band we are achieving 
~20 mmag and ^--^30 mmag RMS accuracy at r ~18 mag for PSF- 
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Figure 4. (a): Top row: Histograms of galaxy proper motions in RA (left) 
and Dec (right) before astrometric recalibration. Bottom row: Density plots 
of galaxy proper motions in RA (left) and Dec (right) as a function of RA 
before astrometric recalibration. The intensity bar has units of number of 
galaxies deg~ ^ mas" ^ yr. (b): The same as Figure [4(a)] but after astrometric 
recalibration. 

like and non-PSF-like objects, respectiveljlfl The RMS diagrams 
for the other wave bands are very similar except for the u band, 
where we achieve ~23 mmag and ~50 mmag RMS accuracy at 
r ~18 mag for PSF-like and non-PSF-like objects, respectively. 

2.4 Further Astrometric Calibrations 

Each SDSS imaging run was astrometrically calibrated against 
the US Naval Observatory CCD Astrograph Catalog (UCAC; 
IZacharias et alj [2OOO), yieldin g absolute posi tions accurate to 
~45 mas RMS per coordinate Ipier et alj|2003l) . The accuracy is 
limited primarily by the accuracy of the UCAC positions (~70 mas 
RMS at the UCAC survey limit of i? ~ 16), as well as the density 

^ The upturn in the RMS amplitude for stars brighter than r --^ 16 is due to 
a combination of factors including: the appearance of detectable asymmet- 
ric low surface-brightness structure (e.g. diffraction spikes) in the images; 
the effect of the large angular size of the images on the determination of the 
sky-background. 
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Figure 5. Plot of the mean offsets in RA and Dec relative to the reference 
run 5823 as a function of RA for runs 94 (solid lines) and 5918 (dashed 
lines). For both runs we show the offsets for camera column 1, and for 
run 94, the lower and upper lines correspond to the RA and Dec offsets, 
respectively. Note that this is but a small section of the full stripe. 



of UCAC sources. The version of UCAC used to calibrate SDSS 
lacked proper motions, thus any proper motions based on SDSS 
positions will be systematically in error by the mean proper motion 
of the UCAC calibrators. 

We illustrate the systematic errors inherent in the SDSS 
astrometry by considering galaxies in the magnitude range 
17 < r < 19.5 that have light-motion curves with at least 
20 astrometric measurements. For this set of galaxies we measure 
the proper motion in RA and Dec, and we find that the galaxies are 
systematically moving with proper motions of the order of 10 mas 
per year in both right ascension and declination! The problem is 
clearly evident in Figure |4(a)| where we show histograms of the 
galaxy proper motions in RA (top left hand panel) and Dec (top 
right hand panel). The trends of galaxy proper motion with RA are 
shown in the bottom left hand panel for proper motion in RA and 
in the bottom right hand panel for proper motion in Dec. 

Proper motions based on multi-epoch SDSS data can be im- 
proved by recalibrating each SDSS imaging run against a reference 
SDSS run, rather than using the UCAC catalogue positions. The ac- 
curacy of th e relative astro metry between runs is ~20 mas RMS per 
coordinate jPier et al ."2003). far superior to the accuracy of the ab- 
solute positions, due both to the more accurate centroids, as well as 
the far greater density of calibrators, for SDSS compared to UCAC. 
Further, by using galaxies as calibrators, the proper motions can be 
tied to an extragalactic reference frame and are thus inertial. This is 
the method used in this paper to correct the systematic SDSS astro- 
metric errors illustrated in Figure |4(a)| However, it is appropriate 
to mention that the positional system of the calibrators still refers 
to the epochs given by UCAC. 

All imaging runs along the north strip have been recalibrated 
using run 5823 as the reference run. All imaging runs along the 
south strip have been recalibrated against run 4203, after first re- 
calibrating run 4203 against run 5823. In order to recalibrate a tar- 
get run, offsets in RA and Dec are calculated for matching "clean" 
galaxies in the reference run (rejecting galaxies affected by prob- 
lems with deblending, pixel interpolation, multiple matches, etc. in 
either run). Only galaxies in the magnitude range 17 < r < 19.5 
are used to avoid large galaxies with poorly defined centroids. For 
each object in the target run, the mean offsets in RA and Dec for the 



nearesjfl 100 such galaxies are calculated and added to the object 
position. This procedure recalibrates the positions in the target run 
to the reference frame defined by the galaxies in the reference run. 

In Figure [5] we show example mean offsets in RA and Dec 
for camera column 1 from runs 94 and 5918. Notice that run 94, 
observed in 1998, requires larger mean offsets to correct for the 
mean proper motion of the UCAC calibration stars than run 5918, 
observed in 2005, since it is further away in time from when the 
reference run 5823 was observed in 2005. 

After recalibrating the astrometry for all light-motion curves 
in the LMCC, we have recreated Figure [4(a)] as Figure [4(b)] using 
the same sample of galaxies. The histograms of the galaxy proper 
motions in RA and Dec are now centred around '^O mas yr~^ in- 
dicating that galaxies are stationary in the recalibrated astrometric 
system of the LMCC. Also, the lower panels demonstrate that the 
RA dependence of the galaxy proper motions has been properly 
removed. 

There is also some evidence that the galaxy proper motion 
scatter has been improved. The RMS deviation of the residuals 
about a fourth-degree polynomial fit in each of the bottom panels of 
Figure [4(a)[ is 5.4 mas yr^^ for RA and 5.2 mas yr~^ for Dec. This 
may be compared to the improved RMS deviation of the proper mo- 
tions in each of the bottom panels of Figure [4(b)[ at 4.8 mas yr~^ 
for RA and 4.6 mas yr~^ for Dec. 

The SDSS pipelines do not supply uncertainties on the mea- 
sured celestial coordinates in the tsObj files, and so we have deter- 
mined a noise model describing how the astrometric noise behaves 
as a function of magnitude. This was done by examining the dis- 
tribution of coordinate RMS for objects in the LMCC. However, 
we found that the astrometric noise in the 2005 observing season 
was noticeably larger than in previous seasons, most likely due to 
the less stringent restrictions on observing conditions leading to a 
greater spread in PSF full-width half-maximum and object signal- 
to-noise. To properly account for this, we determined separate noise 
models for the pre-2005 and 2005 observing seasons. 

To determine the astrometric noise models, we select all PSF- 
like objects (MEAN.OB JECT.TYPE = 6) with at least 20 good 
epochs in r in each of the pre-2005 and 2005 observing seasons. 
For these objects we derive the distribution of coordinate RMS de- 
viations for 0.5 mag bins for both pre-2005 and 2005 data, and 
fit a peak and dispersion for each bin. In Figure |6(a)[ we plot the 
peak RA RMS deviation for each magnitude bin versus r mag- 
nitude for pre-2005 data (filled circles) and 2005 data (open cir- 
cles; offset by 0.15 mag to the left for clarity). We obtain very 
similar results for the Dec coordinate. We fit the peak data as a 
function of magnitude m via an exponential function of the form 
/(m) ^ A + B exp (C (m - 18)) where ^, B and C are fitted 
parameters, and plot the fitted models in Figure [6(a)] as continuous 
and dashed curves for pre-2005 and 2005 data, respectively. 

The following equations represent our final adopted astromet- 
ric noise model, based on the exponential model fits for both the 
RA and Dec coordinates: 

Pre-2005: (J<,(t) = as{t) = 32.0 + 0.430 exp (1.34 {m{t) - 18)) 
2005: ac.it) = as{t) = 35.4 + 0.783 exp (1.09 (m(t) - 18)) 

(1) 

where (Ja{t) and as{t) are the uncertainties on the measured ce- 

Nearest in coordinate parallel to the scan dii'ection, right ascension - the 
coordinate perpendicular to the scan direction, declination, is ignored, as 
the length of the binning window is always larger than the width of a scan. 
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Figure 6. (a): Plot of the peak RA RMS deviation in 0.5 mag bins versus r magnitude for PSF-like objects (stars) for pre-2005 data (filled circles) and 
2005 data (open circles; offset by 0.15 mag to the left for clarity), (b): Plot of the peak RA RMS deviation in 0.5 mag bins versus r magnitude for non- 
PSF-lrke objects (galaxies) for pre-2005 data (filled circles) and 2005 data (open circles; offset by 0.15 mag to the left for clarity). Both panels: Error bars 
represent the dispersion in the distribution of RMS deviations in each magnitude bin. Plotted functions (continuous and dashed curves) are all of the form 
/(m) = A + B exp (C (m — 18)) where m denotes magnitude and A, B and C are fitted parameters. 



lestial coordinates a{t) and 5{t), respectively, at time t, and m{t) 
represents the brightest PSF magnitude out of the five photomet- 
ric measurements at time t. Evidence that this noise model is valid 
comes from the fact that the distribution of per degree of free- 
dom of the proper motion fit for the HLC (Section 3. 1) is peaked at 
a value of ~ 1 . 1 . Note that astrometric uncertainties are not given in 
the LMCC and should be obtained via Equation [T] 

In Figure |6(b)[ we plot the results of the same co- 
ordinate RMS deviation analysis for non-PSF-like objects 
(MEAN_OBJECT_TYPE = 3) with at least 20 good epochs in r 
in each of the pre-2005 and 2005 observing seasons. It is clear 
from the plots in Figure |6] that we are achieving ~32 mas and 
~35 mas RMS accuracy at r ~18 mag for stars for pre-2005 and 
2005 data, respectively, and ~35 mas and ~46 mas RMS accuracy 
at r ~18 mag for galaxies for pre-2005 and 2005 data, respectively. 

2.5 Catalogue Format 

The LMCC exists as eight tar files, one for each hour in 
RA from 20*^ to 4*'. Each tar file contains 60 subdirecto- 
ries corresponding to the minutes of RA, and the light-motion 
curves are stored in these directories based on their mean RA 
coordinates. The LMCC contains 3700548 light-motion curves, 
2807047 of which have at least 20 epochs. The tar files 
(~29.5 Gb compressed) may be obtained by web download from 
http://das.sdss.org/value_added/stripe_82_variability/SDSS-82.public/ 
Light-motion curve plotting tools written in IDL may also be 
downloaded from the same website. 

A single light-motion curve is stored as an ASCI I file with a 
name constructed from the unweighted mean position of the cor- 



responding object. The ASCII light-motion curve file contains a 
header line describing the column meanings, followed by exactly 
five rows for each epoch (one row for each wave band) in strict 
time order. All five wave band measurements are included for com- 
pleteness, even though it is possible that at any one epoch, up to 
four wave band measurements may not satisfy the quality criteria 
described in Section 2.3. In Table[3]we describe the columns that 
make up a light-motion curve from the LMCC. 

Figure [7] shows some clear examples of photometric vari- 
ability and motion from the LMCC. Figure |7(a)| presents the 
lightcurve in r (upper points) and g (lower points) of the large- 
amplitude long-period variable star SPSS J2205 14 .58-1-000845.7, 
most likely a Mira variable ( Watk ins et alj |2008|) . Figure |7(b)| 
presents the motion curve of the known ultr acool white dwarf 
SDSS J224206. 19-1-004822.7 iKilic et al.ll200^ . Both panels illus- 
trate the dramatic increase in temporal sampling produced by the 
start of the SDSS-II Supernova Survey in 2005. 



3 THE HIGHER-LEVEL CATALOGUE 
3.1 Catalogue Description 

The HLC supplies a set of 229 derived quantities for each light- 
motion curve in the LMCC. These quantities are aimed at describ- 
ing the mean magnitudes, photometric variability and astrometric 
motion of the objects in the LMCC, and they are calculated using 
only light-motion curve entries that satisfy the quality constraints 
from Section 2.3. Those quantities in the HLC related to photom- 
etry are described in Table [4] while those related to astrometry are 
described in Table[5] 
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Table 3. The list of columns that make up a light-motion curve in the LMCC along with a brief description. 



Column Nuinber Coluinn Naine 



Type 



Description 



1 


Run 


INTEGER 


SDSS imaging run 


2 


Rerun 


INTEGER 


Version of the SDSS frames pipeUne used to process the data 


3 


Field 


INTEGER 


SDSS field number along a strip 


4 


Camcol 


INTEGER 


SDSS camera column 


5 


Filter 


INTEGER 


SDSS wave band (0 = u, 1 = g, 2 = r, 3 = i, 4 = z) 


6 


Object Type 


INTEGER 


Object classification (3 = Galaxy, 6 = Star) 


7 


RA 


DOUBLE 


Right ascension J2000 (deg) 


8 


Dec 


DOUBLE 


Declination J2000 (deg) 


9 


Row 


FLOAT 


CCD row coordinate (pix) 


10 


Coluirm 


FLOAT 


CCD column coordinate (pix) 


11 


HJD 


DOUBLE 


Heliocentric JuHan date (days) 


12 


PSF Luptitude 


FLOAT 


PSF magnitude (lup)* 


13 


PSF Luptitude En'or 


FLOAT 


Uncertainty on the PSF inagnitude (lup)* 


14 


PSF Flux 


FLOAT 


PSF flux normalised by the flux from a zero-th magnitude object 


15 


PSF Flux Error 


FLOAT 


Uncertainty on the PSF flux 


16 


Exp Luptitude 


FLOAT 


Exponential magnitude (lup)* 


17 


Exp Luptitude Error 


FLOAT 


Uncertainty on the exponential magnitude (lup)* 


18 


Exp Flux 


FLOAT 


Exponential flux normalised by the flux from a zero-th magnitude object 


19 


Exp Flux Error 


FLOAT 


Uncertainty on the exponential flux 


20 


Sky 


FLOAT 


Sky background brightness (lup arcsec"^)* 


21 


Sky EiTor 


FLOAT 


Uncertainty on the sky background brightness (lup arcsec"-^)* 


22 


FWHM 


FLOAT 


Full-width half-maximum of the PSF (arcsec) 


23 


Photometric Calibration Tag 


INTEGER 


Flag indicating the photometric calibration status (1 = Calibrated, = Uncalibrated) 


24 


Photometric Zero Point 


FLOAT 


Fractional flux offset applied to the flux values 


25 


Flag 1 


LONG 


Object flags (tsObj file tag FLAGS) 


26 


Flag 2 


LONG 


More object flags (tsObj file tag FLAGS2) 


27 


Astrometric Calibration Tag 


INTEGER 


Flag indicating the astrometric calibration status (1 = Calibrated, = Uncalibrated) 


28 


RA Correction 


DOUBLE 


Correction applied to right ascension (deg) 


29 


Dec Correction 


DOUBLE 


Correction applied to declination (deg) 



"See lLupton. Gunn & Szalavl h999l) and lStoughton et alj j2002l) for flie definition of a luptitude. 



In Table |4l if a tag name is associated with a 5-element ar- 
ray, then the 5 values represent the described quantity for each of 
the five SDSS wave bands in the order u, g, r, i and z. When a 
certain wave band has insufficent "good" light-motion curve en- 
tries to calculate a particular quantity, a value of zero is stored (this 
also applies to Table [Sjl. For instance, the first value in the array 
MEAN_PSFMAG is set to zero for any light-motion curves with no 
"good" entries for the u band. 

All quantities in Table|4]with CLIP at the end of the tag name 
are calculated using a 4o--clip algorithm that rejects only the worst 
outlier at any one iteration, and terminates when no more outliers 
are identified. Similarly, all quantities in Table|4]with ITER at the 
end of the t ag name are ca lculated using the iterative procedure 
described in IStetsonI ( Il996 l) to dynamically reweight data points 
based on the size of the residuals from the mean. Both these sets 
of quantities have been designed to be more robust against outliers 
than a simple inverse variance weighted mean. 

The SDSS photometric pipeline performs a morphological 
star/galaxy separation, the quality of which is intimately related 
to seeing and sky brightness. While the accuracy is very good for 
bright objects, there can be confusion for faint objects. The quan- 
tity MEAN.OB JECT.TYPE in Table|4]is an unweighted mean of the 
SDSS object type classification. Hence it has a value of 3 if the ob- 
ject is classified as a galaxy at all epochs, a value of 6 if the object 
is classified as a star at all epochs, and a value between 3 and 6 
otherwise. The reliability of MEAN.OB JECT.TYPE for object type 
classification depends on the reliability of the SDSS object type 
classifier and the number of epochs at which the object was ob- 
served. 



The quantity MEAN.CHILD in Table|4]is an unweighted mean 
of whether the child object bit is set or not. In other words, this 
quantity has a value of 1 if, at all epochs, the object results from the 
deblending of a parent object, a value of if the object was never 
the result of the deblending of a parent object, and a value between 
and 1 otherwise. 

The Galactic extinctions described by EXTINCTION in Ta- 
ble [4] are derived us ing the maps of dust column density from 
iSchlegeletalJ l ll998h . 

For eclipse and flare detection in light curves we include the 
statistics ECL.STAT and FLARE.STAT (Table|4j in the HLC. In 
calculating these values, we assume that any eclipse/flare event in 
a light curve will include only one photometric data point since 
the time elapsed between consecutive scans is at least one day. 
Hence, for each photometric data point (using PSF magnitudes 
only), we calculate the statistic S , based on a matched filter from 
iBramich et al.. (2005,) . and defined by: 



Xconst 



Xout 



(Xout/z^) 



(2) 



where Xconst chi-squared value of a constant fit for each wave 
band to the whole light curve, and Xout '^^ the chi-squared 

and number of degrees freedom, respectively, of a constant fit for 
each wave band to the out-of-eclipse/flare light curve. In order to 
avoid false positives, we only calculate for epochs with photo- 
metric data points that are "good" in at least two wave bands. The 
adopted values of ECL.STAT and FLARE.STAT are then taken 
to be the largest values of 5*^ for photometric data points fainter 
and brighter, respectively, than the mean. We also record the corre- 
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Table 4. The list of derived quantities related to photometry that are stored for each light-motion curve in the HLC. These quantities are calculated using only 
light-motion curve entries that satisfy the quality constraints from Section 2.3. 



Tag Name In HLC 


Type 


Description 




STRING 


T iCTnt-mofinn r'lin/f filfnfimp 




STRTNn 


WUjCLl IlcllllC 111 iJi_-/iJ>j LJaVa rvClCaftC U ^lllLClllclliiJllill rt.ftlHJlHJllllLcll LJ lllUll tllJlJUJVCU. HJillld.1^ 




5 V TNTFGFR 


^Tiimhpr ni ar\nc\ nhritnmptnr" Hntn nr\intQ 

i 1 LIlllU^l yjx ^\J\J\A yjLl\J\\jLLi\^ 11 IV^ UuLu IJVJllJLj 




5 V FT OAT 


TnvprQP vnfianr'P wpiCTntpn mprin ot tVip T'SF mntrnitnHpc 

ltIVCl3& Vullull^^ Wd^IIL^Ll JllCClil \J1 LiiC 1 l31 illu^Ill LLiLl^d 


MF7\M pi^FMar FRR 


S V FT OAT 


T Tnr-prtiiinfv nn MFAM P^iFMAP 


MEAN EXPMAG 


5 X FLOAT 


Inverse variance weighted mean of the exponential magnitudes 




S V FT OAT 


TTnrprtnintv nn MFAM FYPMAP 
u iiuci LdiiiLy uii 1 iiLirtiM _rj A jr i Jrto 




5 X FT OAT 

A 1. XjV>/^1 


T? r»nt-mpnn-<;niinT"P Hpvintinn r\T thp T^SF mnonitiinpc 

IV^J^L IJlCClll^U LluL ^ \Xk^ V lu.LlWiI VJl LiiC 1 Ol illu^Ill LLiLi^A 


RM^^ pypMAr; 


5 V FT OAT 


T? r»nl"-mpnti-Qniinrp Hpvintinti r\\ tVip PYnnnpntial mno"nihiHpQ 

IVU^JL IJlCtlll^U LLu.1 X-iC V lu.Li\JiI VJl LliC ^.j\L/\Jll^llll£ll 111£1^1J1LL1\J.^A 




S V FT OAT 


l^ni _cni 1 nrpri r\T tnp T'^F mifmitnripc 

aULlCliCU Ul lllC 1 vJl lIldgllllULIC^ 


THmn FypMA(^ 


5 X FT OAT 


f^ni -QfiiiarpH nf thp PYnnnpntinl mnCTmtiirip*; 

V^IIl^VJ LIul ^VJ. ^1 LlJ^ &.j\LJ^il^iILiul Jlia^lllLLlLl&d 


M nnnn FPriPHc; pqt? tt tp 


S V TNTFHFT? 


i> uiiiuci \ji gijuii r oi iiicigiiiiutics di ici 'tc v^iiuuiiig 


IN _oUUJJ_Ej IT no —111 A Jr ij ± r 


S V TNTFnFT? 


Number of good exponential magnitudes after ^tr -clipping* 




^ V FT OAT 


^jT«r*l inr^pH in\7Prcp \/^iti iinr'P \x7PiO"htpH mp^in £\x tHp T-^^F rnQfrnitiiriPC^ 
'tC UllUUCU lllVCloC VullollL.C' WClgHlCU lllCall \jl LllC JTOl lllaXII^^LULlCd 


MFAW PtUFMArt FRR TT TP 


S V FT OAT 


TTnrprtninh? r»n MFATvT pciFMAP PT TP* 
u iiv^ci iciiiiiy uii r'iHiri.j.N tr or i*irio_i.^Xj x jr 


MFAKi FypMAr; CT TP 


5 X FT OAT 

^ A J. I ^\ J r \ 1 


^iT-plinnp/H in^'PFQP vnrinnpp wpioHtpH mpnn nf thp pvnnnpntinl mricTnitiirlpc* 

^(-/^lllJlJCVj lllVdL^t^ V Ol lall^C W t^l^llLCU lll^ail VJl LllC CA.IJVJlldlLlu.1 111£1^111LLI1-ICl) 


MEAN_EXPMAG_ERR_CL I P 


5 X FLOAT 


Uncertainty on MEAN_EXPMAG_CLIP* 


■RM'^ pt^FMAP PT TP 
rs.i 1 o _jr o r i Jrt.o_v. j_j ± jr 


^ V FT OAT 


T< oj^l"_m pnn«cn 1 1 ni"p rlp\^iif'ir^n c^x trip ^i^_r*Mr\npri T'nIF ■mcmnitiiripc 

IVUtJL lllCtlll aUUttlC UCVltlLltJll \Jl LllC 'T^fJ CllUUCLl ITiJi llldglllLLILlCa 


RMC; FYPMAP PT TP 


^ V FT DAT 


t< o/^l"_m pcin«cn 1 1 qtp ri p\^i nti rMi c^x tn p li trnpri pvnr^npn f"icil m ci trn iti iHpc 

IVUUL lllCctll SUUttlC UCVltlLltJll \JL LllC CllUUCU CAUtJllCllLldl lllaglllLUUCft 


ppjTC^P) P'^FMAr; fT.TP 


5 X FT OAT 


r^ni -'iniifn'prl r\\ thp ^iT-plinnprl T'SF mfionituHp*;* 

V^lll L>U LItll Ctl VJl LllC ClllJIJCU 1 ijl lllu^lllLLILlC3 


fUTcn FypMAr; fT.TP 


5 X FT OAT 

^ A J. I ^\ J r \ 1 


r^ni -';niifii'p(H r\\ thp 4/T-flinnprl pvnnnpntiril maonitiirip';* 

V^lll L^U LIul C^J yji. LllC ClllJIJCLi CA.IJVJllCllLltll 111£1^111LLIU.Cl> 


MFAM PCIFMAP TTFR 


^ V FT OAT 


TtpriifpH iTT\^PTCP ■\/'iin ^inr*p as/pi trhtpfi mpcin c^x thp T^xF mntmitiiripc 
ILClaLCLl lllVCl&C VallcUICC WClglllCU lllCall \jl LllC JT l31 lllu^lllLLlLICA 


MFAN P^FMAr; FRR TTFR 


5 X FT OAT 


TTnrprtaintv nn MFAN P'^FMAr; TTFR* 

\^ llCCl LdlllLy ^Jll I"iJ_jri.i\ _r o r i Irto— X X J_j Jr\. 


M'PZiW PYPMSP TTPR 
FiJljiiiN_rj Airi\LH.o L 1 i!jJ:\ 


S V FT OAT 


Iterated inverse variance weighted mean of the exponential magnitudes 


MFAKI FYPMAn FRR TTFR 


S X FT OAT 


TTnrprtaintv nn MFAN FYPMAi^ TTFR* 

llV^Cl Lu.lllLy KJLi 1 IHirtlN _£_] A jr L J_rt.O_ X XI_jX\ 


Pp-DPpM'TTT.p nS PC:F 


S V FT OAT 


Sth T'pri^pntilp OT thp r*iimiiliti\/p /iiclrirMitmn or T'vF mnnnitiiripc 

JLll JT Cll-ClllllC Ui LllC ^LIlllLlluLl VC Uldll lUULlUll Ui IT i31 lllaglllLULlC A 


p-frpp-piyTTTTF Rfi Pc;F 


^ V FT OAT 


A/jpriiin r\T thp T-'vF mfitrnitiiripc 

iVlCUltlll UI LllC JTiJl Illd^lllLLIVlCA 


PFRPFMTTTF pc;F 
irHiiAL-iijiM X xxiiL J -J _ir o r 


^ V FT DAT 


QSth T'pT'r^pn ti 1 p c^x thp r'iimiilQti\^p riictnhiitmn /tF T-'vF ■mQcrnitii/Hpc 
7 J 111 IT CI L.C11111C \JL LllC L.LIlllLllaLlVC Ula LI lULlLlUll UI JTOl lllaglllLLUJCa 


pppr-pia-TTT -f? ns FYP 
IT xlj r\ ^iZji^i X xxjiIj_u ^ /i. Jr 


^ V FT DAT 


Sth T-'pf/^pn flip OT thp r*i imi 1 1 nti^^p rii ctn hi iti nn nT pYt^nnpntinl m ntrnituHpc 

JLll JT ClL/CllLllC Ui LllC UUIIILllaLI V C Ulft LI lULlllUll UI CA UUIlCllLltll llluXII^^ LL1L1C& 


PFPPFKTTTT F Rfl FYP 

xILrxL-IljiN X XXjI1j_J W tlj A. IT 


^ V FT OAT 


A/Iprii^in nt thp PYT\nnpnticil TncitmituHpc 

iVlCLllClll Ui LllC CAL/UllCllLltll IIIClglllLLILlCd 


PFRPFMTTTF QR TTYP 
iriijrs.v^iijiN X xxjI1j_!3 j tij a Jr 


S V FT OAT 


95th Percentile of the cumulative distribution of exponential magnitudes 


T T MF c; P A M 
X X i^IJZj _ O x rtiN 


FLOAT 


T^imp crxiin nt thp h frht-mnti nn r*i in?p i H » 

1 illlC aUtlll Ul IIIC liglll IIIULIUIl CLII VC \ 




FLOAT 


Unweighted mean of the object classification 


MFAN rHTT.n 
i*XLjr^Nj n X xju 


FLOAT 


T Tn\x/pi chtpH mpnn nf whpthpr thp philH hit ic QPt ni" nnt* 

11 W Cl^llLCvl lllCull Ul WllCLIlCl LllC UllllLl UlL Ic^ jCL Ul IIUL 


TT'YTTMPTTnM 


S V FT OAT 


Galactic extinction (mag) 


ECL.REDCHISQ.OUT 


FLOAT 


Reduced chi-squared out-of-eclipse for the PSF magnitudes* 


ECL.STAT 


FLOAT 


Eclipse statistic for the PSF magnitudes* 


ECL.EPOCH 


DOUBLE 


Eclipse epoch as a heliocentric Julian date (d)* 


FLARE .REDCHISQ.OUT 


FLOAT 


Reduced chi-squai'ed out-of-flai"e for the PSF magnitudes* 


FLARE.STAT 


FLOAT 


Flare statistic for the PSF magnitudes* 


FLARE .EPOCH 


DOUBLE 


Flare epoch as a heliocentric Julian date (d)* 


STETSON.INDEX.J.PSF 


5 X FLOAT 


Stetson J-index for the PSF magnitudes* 


STETSON.INDEX.J.EXP 


5 X FLOAT 


Stetson J-index for the exponential magnitudes* 


STETSON.INDEX.K.PSF 


5 X FLOAT 


Stetson K-index for the PSF magnitudes* 


STETSON.INDEX.K.EXP 


5 X FLOAT 


Stetson K-index for the exponential magnitudes* 


STETSON.INDEX.L.PSF 


5 X FLOAT 


Stetson L-index for the PSF magnitudes* 


STETSON.INDEX.L.EXP 


5 X FLOAT 


Stetson L-index for the exponential magnitudes* 


VIDRIH.INDEX.PSF 


5 X FLOAT 


Vidrih index for the PSF magnitudes* 


VIDRIH.INDEX.EXP 


5 X FLOAT 


Vidrih index for the exponential magnitudes* 


*See text for more detail. 
tjVIay be empty. 



spending values of Xout/'^ f°r ECL.STAT and FLARE.STAT in the 
quantities ECL.REDCHISQ_OUT and FLARE_REDCHISQ_OUT re- 
spectively, along with the epoch of the putative eclipse/flare event 
in ECL.EPOCH and FLARE_EPOCH respectively. 

^In Table [4] the Stetson variability indices J, K and L jStetsorJ 

Il996h for both the PSF and exponential magnitudes are stored in 
the quantities with tag names starting STETSON.INDEX. We chose 
the r band as the comparison wave band and consequently the third 
element in each of these six quantity arrays is set to zero. 



Our final measure of light curve variability is via a quan- 
tity called the Vidrih variability index. RMS magnitude devia- 
tions for non-variable lightcurves plotted versus mean magnitude 
m in a given wave band (referred to as an RMS diagram) are 
scattered around a three parameter empirical exponential func- 
tion /(m) = A + B exp (C (m - 18)), while 
the RMS magnitude deviation of variable sources is expected to 
be noticeably larger. We consider two lightcurve samples, that of 
PSF-like objects (MEAN.OBJECT.TYPE = 6) and non-PSF-like ob- 
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Table 5. The list of derived quantities related to astrometry that are stored for each light-motion curve in the HLC. These quantities are calculated using only 
light-motion curve entries that satisfy the quality constraints from Section 2.3. 



Tag Name In HLC 


Type 


Description 


RAJIEAN 


DOUBLE 


Inverse variance weighted mean of the RA measurements (deg) 


RAJ4EAN.ERR 


FLOAT 


Uncertainty on RAJ4EAN (deg) 


RA.PM 


FLOAT 


Proper motion in the RA coordinate (arcsec yr~^) 


RA_PM_ERR 


FLOAT 


Uncertainty on RA.PM (arcsec yr~^) 


RA.CHISQ.CON 


FLOAT 


Chi-squared of the RA measurements for a model that includes only a mean position 


RA.CHISQ.LIN 


FLOAT 


Chi-squared of the RA measurements for a model that includes a mean position and a proper motion 


RA_MEAN.CLIP 


DOUBLE 


Inverse variance weighted mean of the clipped RA measurements (deg) 


RA J4E AN.ERR.C L I P 


FLOAT 


Uncertainty on RAJ4EAN.CLIP (deg) 


RA.PM.CLIP 


FLOAT 


Proper motion in the RA coordinate for the clipped RA measurements (arcsec yr~^) 


ra.pmj:rr.clip 


FLOAT 


Uncertainty on RA.PM.CLIP (arcsec yr""*^) 


RA.CHISQ.CON.CLIP 


FLOAT 


Chi-squared of the clipped RA measurements for a model that includes only a mean position 


RA.CHISQ.LIN.CLIP 


FLOAT 


Chi-squared of the clipped RA measurements for a model that includes a mean position and a proper motion 


DECJ^IEAN 


DOUBLE 


Inverse variance weighted mean of the Dec measurements (deg) 


DECJIEAN.ERR 


FLOAT 


Uncertainty on DECJ^EAN (deg) 


DEC.PM 


FLOAT 


Proper motion in the Dec coordinate (arcsec yr~ ^ ) 


DECJ'M.ERR 


FLOAT 


Uncertainty on DECJ'M (arcsec yr^^) 


DEC.CHISQ.CON 


FLOAT 


Chi-squared of the Dec measurements for a model that includes only a mean position 


DEC.CHISQ.LIN 


FLOAT 


Chi-squared of the Dec measurements for a model that includes a mean position and a proper motion 


DECJXIEAN.CLIP 


DOUBLE 


Inverse variance weighted mean of the clipped Dec measurements (deg) 


DEC31EAN.ERR.CLIP 


FLOAT 


Uncertainty on DECJ^EAN.CLIP (deg) 


DECJ'M.CLIP 


FLOAT 


Proper motion in the Dec coordinate for the clipped Dec measurements (arcsec yr~ ^ ) 


DEC_PM.ERR.CLIP 


FLOAT 


Uncertainty on DECJ'M.CLIP (arcsec yr~^) 


DEC.CHISQ.CON.CLIP 


FLOAT 


Chi-squared of the clipped Dec measurements for a model that includes only a mean position 


DEC.CHISQ.LIN.CLIP 


FLOAT 


Chi-squared of the clipped Dec measurements for a model that includes a mean position and a proper motion 


TO 


DOUBLE 


Inverse variance weighted mean of the heliocentric Julian dates using the uncertainties on the astrometric measurements (d) 


TO.CLIP 


DOUBLE 


Inverse variance weighted mean of the heliocentric Julian dates for the clipped astrometric measurements (d) 


N.POS.EPOCHS 


INTEGER 


Number of good positional measurements 


N.POS.EPOCHS.CLIP 


INTEGER 


Number of good positional measurements after clipping 



jects (MEAN.OBJECT.TYPE = 3). For each lightcurve sample and 
wave band, we iteratively fitted /(m) to the corresponding RMS 
diagrams using a 3(T-chp algorithm, employing PSF magnitudes 
for the star lightcurve sample and exponential magnitudes for the 
galaxy lightcurve sample. In each case we also constructed a func- 
tion g{m) describing the standard deviation of the scatter around 
/(m) via a four-degree polynomial fit to the standard deviation of 
the RMS magnitude deviations measured in 0.25 magnitude bins. 
Then, for any light curve, the Vidrih index V was calculated as its 
RMS magnitude deviation minus f{m), normalized by g{m), with 
negative values set to zero, and it is stored in the quantities with 
tag names starting VIDRIH.INDEX (Table|4ll. Note that due to the 
way the Vidrih indices are constructed, VIDRIH.INDEX_PSF is 
only relevant to stars and VIDRIH. INDEX_EXP is only relevant to 
galaxies. 

In Table [5] we describe the quantities associated with the ob- 
ject astrometry. Specifically we fit a proper motion model for ce- 
lestial coordinates {a, 5), in degrees, at heliocentric Julian date t, 
in days, for each light-motion curve: 

7.605 X 10~^RA.PM , 

a(t) = RAJ1EAN H ^ — (t - TO) (3) 

^ ' cos(DEC_MEAN) ^ ' 

S{t) = DECJ1EAN + 7.605 x lO'^DECJM {t - TO) (4) 
We set TO as the weighted mean of the epochs of observation: 



TO 



E,i/'^cfe)2 



(5) 



where ti represents the set of heliocentric Julian dates for the posi- 
tional measurements, and aa{ti) represents the set of uncertainties 



on the measured celestial coordinates a,; and Si . The uncertainties 
aa{ti) are calculated via Equation[T] 

We solve for the quantities RAJIEAN, RA.PM, DEC31EAN 
and DECPM by minimising the appropriate chi-squared using a 
downhill simplex algorithm. We calculate the uncertainties on these 
quantities, namely RA.MEAN_ERR, RA_PM.ERR, DECJ4EAN.ERR 
and DEC_PM_ERR, by assuming that cos(DEC.MEAN) 1 (a 
valid assumption for \S\ < 1?26) and noting that minimis- 
ing the chi-squared for each of the resulting equations is a lin- 
ear least-squares problem with an analytic solution. We supply 
the best-fit chi-squareds via the quantities RA_CHISQ_LIN and 
DEC.CHISQ.LIN, and we record the number of positional mea- 
surements used in the fit in the quantity N.POS.EPOCHS. We also 
supply the chi-squareds RA.CHISQ.CON and DEC.CHISQ.CON 
for a model including only a mean position, which facilitates the 
calculation of the following delta chi-squared: 



Axl = RA.CHISQ.CDN - RA.CHISQ.LIN 

+ DEC.CHISQ.CON - DEC.CHISQ.LIN 



(6) 



The statistic Axfi may be used to calculate the significance of a 
proper motion measurement by noting that it follows a chi-square 
distribution with two degrees of freedom. 

The whole fitting process is also iterated using a clipping al- 
gorithm that rejects the worst outlier at any one iteration, and ter- 
minates when the change in the fitted proper motion is smaller than 
3 mas yr~^. Consequently, all the derived astrometric quantities 
described so far have a corresponding quantity calculated for the 
clipped positional measurements and stored in the HLC (all quanti- 
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Figure 7. (a): Light curve of the Mira variable star candidate 
SDSS J2205 14.58+000845.7 in the r (upper points) and g (lower points) 
bands, (b): Fitted proper motion in RA (track at bottom left, left-hand y- 
axis) and Dec (track at top left, rig ht-hand y-a xis) for the ultracool white 
dwarf SDSS J224206. 19+004822.7 tolic et al.ll2006,) . 

ties listed in Table[5]with a tag name ending CLIP). The astromet- 
ric quantities derived from the chpped astrometric data are more 
robust than those derived from the undipped astrometric data. We 
strongly advise that any statistical studies undertaken with the as- 
trometry in the HLC should only use the clipped quantities. 

Finally, we have attempted to detect the parallax signal for 
nearby stars by fitting a parallax model to each light-motion curve. 
However, our investigation has made it clear that the distribution 
of observations during the same few months each year along with 
the positional accuracy of the astrometric measurements are not 
enough to enable the detection of a clean parallax signal. 

3.2 Catalogue Format 

The HLC is stored in eight FITS binary tables, one for each hour 
in RA from 20*' to ^ . Each record in the HLC corresponds to a 
single light-motion curve and stores 229 derived quantities with 
tag names listed and described in Tables |4] and |5] The FITS files 
(~1.9 Gb compressed) may be obtained by web download from 
^http://das.sdss.org/value.added/stripe^82j^>ariability/SDSS-82.public/ 

The IDL Astronomy User's Librarjlf] function "mrdfits" is a 
convenient way to read-in the HLC, storing the FITS binary table 

^ The IDL Astronomy User's Library is web-hosted at 
|http://idlas tro.gsfc.nasa.gov/ and maintained by Wayne Landsman at 
the Goddard Space Flight Centre. 



automatically in an IDL structure. The IDL "where" function then 
becomes a very powerful tool for accessing subsets of data with 
ease. 



4 EXTERNAL COMPARISONS 
4.1 Photometry 

The most relevant comparison of our cat alogue photometry can be 
made by comparing our results to those of llvezic et al.l ( 120071) (from 
now on IV07) who use 58 pre-2005 SDSS-I imaging runs of Stripe 
82, observed under mostly photometric conditions, to construct a 
standard star catalogue. IV07 match unsaturated point sources sat- 
isfying high signal-to-noise criteria (photometric uncertainties be- 
low 0.05 mag) between runs and choose sources with at least 4 
epochs. Those objects that are identified as non-varying, by requir- 
ing a x"^ per degree of freedom for the mean magnitudes of less than 
3 in the g, r and i bands, are chosen as candidate standard stars. 
The standard star catalogue is then internally recalibrated, correct- 
ing for flat-field errors and time-variable extinction (for example, 
due to clouds). 

The main differences between the LMCC and the IV07 cata- 
logue are that we use a different photometric recalibration method, 
but based on the same idea of spatially dependent photometric ze- 
ropoints, and that we include the generally non-photometric obser- 
vations from the SDSS-II supernova runs. Not surprisingly, for the 
886396 objects from the IV07 catalogue in the overlap area with 
the LMCC, we find unambiguous matches in the LMCC for 886191 
objects using a match radius of 0.7", where the missing fraction of 
~0.02% can be explained by the different criteria used to construct 
the two catalogues. Of the 886191 matching objects in the LMCC, 
219627, 606658, 879766, 886028 and 856456 objects have at least 
3 good epochs in the u, g, r, i and z bands, respectively, and mean 
PSF magnitudes brighter than 21.5, 21.5, 21.5, 21.5 and 20.5 mag, 
respectively. It is these objects that we use when constructing the 
histograms in Figure [8] 

In Figure [Heach row of panels corresponds to a different wave 
band and the order employed is u, g, r, i and z from the top row 
to the bottom row. The left hand column of plots are normalised 
histograms of per degree of freedom for the mean PSF mag- 
nitudes in the HLC, and they serve to confirm that the photomet- 
ric uncertainties in the LMCC (inherited from the SDSS pipelines, 
with some adjustment during photometric recalibration) are cor- 
rect, since the histograms peak at a value of / {N^ ~ 1) 1- 

The middle column of panels in Figure [8] are normalised his- 
tograms of the uncertainties on the mean PSF magnitudes in the 
HLC (solid lines) and on the mean magnitudes from the IV07 cata- 
logue (dashed lines). It is clear that the mean magnitudes quoted in 
the HLC are approximately twice as precise as those from IV07, 
which is due to the greater number of epochs included in the 
LMCC. However, systematic errors are still likely to dominate the 
mean magnitudes at the ~1% level, which corresponds to the level 
of the photometric recalibrations and to the systematic error intro- 
duced by the slightly different band passes of each column of detec- 
tors in the SDSS camera. We do not correct the LMCC magnitudes 
for the different bandpasses in contrast to IV07. 

The right hand column of panels in Figure [8] are normalised 
histograms of the difference between the HLC mean PSF magni- 
tudes and the IV07 catalogue mean magnitudes. A Gaussian has 
been fitted to each histogram (dotted lines), and the fitted cen- 
tre and sigma are quoted in each panel. The histograms are cen- 
tred around zero at the millimagnitude level, except for the u band 
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Figure 8. Each row of panels corresponds to a different wave band and the order employed is u, g, r, i and z from the top row to the bottom row. The LMCC 
objects used in the plots are matches with the IV07 standard star catalogue and have at least 3 good epochs with mean PSF magnitudes brighter than 21.5 mag 
in the u, g, r and i bands, and brighter than 20.5 mag in the z band. Left hand column: Normalised histograms of per degree of freedom for the mean PSF 
magnitudes in the HLC. Middle colunm: Normalised histograms of the uncertainties on the mean PSF magnitudes in the HLC (soHd lines) and on the mean 
magnitudes from the rV07 catalogue (dashed Unes). Right hand colunm: Normalised histograms of the difference between the HLC mean PSF magnitudes 
and the rV07 catalogue mean magnitudes, along with a fitted Gaussian (dotted hue). 
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Figure 9. Plots of the running 3(j-clipped mean diiference between the HLC proper motions and those of GK04 (middle points) as a function of RA using a 
1° window (upper panels), as a function of Dec using a 0?01 window (iniddle panels), and as a function of mean r PSF magnitude using a 0.1 mag window 
(lower panels). The left hand panels con'espond to proper inotion in RA and the right hand panels correspond to proper motion in Dec. The upper and lower 
sets of smaller points in each panel represent the mean difference plus or minus the running 3cr-clipped standard deviation, respectively. 



where the IV07 photometry is slightly offset from the HLC pho- 
tometry by ~4 mmag. The histograms have sigmas of 28, 13, 11, 
10 and 16 mmag in u, g, r, i and z, respectively, which are consis- 
tent with the scatter in the photometric zeropoints of the standard 
SDSS runs derived in Section 2.2, and with the scatter in the similar 
internal photometric recalibrations of IV07. 

The data set used to construct the IV07 standard star cata- 
logue is a subset of the data used to construct the LMCC, and 
therefore we can use our superior temporal coverage to identify 
standard star candidates that are actually photometrically variable. 



IWatkins et al.l ( l2008h use the HLC to systematically identify vari- 
able stars in Stripe 82, and we choose to follow all but one of their 
cuts on the HLC quantities. Of the 886191 objects from the IV07 
catalogue that have matching objects in the LMCC, 878172 have 
at least 1 1 good epochs in the g band, allowing for the calculation 
of re liable variability and object type indicators. The requirement 
from IWatkins et all ( |2008|) that an LMCC object has a mean ob- 
ject type greater than 5.5 ensures that chosen variables are PSF-like 
(stars). However, this implies that 16529 objects are probably non- 
PSF-like, corresponding to a fraction of ~1.9%. Hence, we do not 
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apply this requirement, and instead we directly apply the remaining 
requirements that a variable object should have a per degree of 
freedom greater than 3 for both the g and r band mean PSF mag- 
nitudes, and a Stetson L index greater than 1 for the g band. We 
find that just 570 IV07 standard stars, from 878172, are variable, 
corresponding to a fraction of ~0.065%. This confirms that even 
with the addition of many more photometric data, the IV07 stan- 
dard stars are in general still found to be constant at the ~0.01 mag 
level. In considering the full IV07 catalogue with 1006849 candi- 
date standard stars, we suspect that ~650 are actually variables. 

4.2 Astrometry 

iGould & KoUmeieJ ( |2004|) (from now on GK04) carefully com- 
bine SDSS Data Release 1 (SDSS DRl) and USNO-B proper mo- 
tions to produce a catalogue of 390476 objects with proper mo- 
tions n > 20 mas yr^^ and magnitudes r < 20 mag. SDSS DRl 

i iroper motions are based on matches of SDSS to USNO-A2.0 
Monet et a i] |l998ll . which suffer from mismatches causing spu- 
rious high proper-motion objects, and systematic trends in proper 
motion. However, by cross-correlating SDSS DRl and USNO-B, 
GK04 successfully removed the vast majority of spurious proper- 
motion stars. Furthermore, by considering the set of spectroscop- 
ically confirmed quasars in SDSS DRl, GK04 calibrate out the 
position-dependent astrometric biases, using a very similar method 
to our recalibration of SDSS astrometry presented in Section 2.4. 

We compare our proper motions in the HLC to those derived 
by GK04 for the 30546 stars from their catalogue with an unam- 
biguous positional match in the HLC using a match radius of 0.7". 
In Figure |9] we plot the running 3(7-clipped mean difference be- 
tween our HLC proper motions and those of GK04 (middle points) 
as a function of RA using a 1° window (upper panels), as a func- 
tion of Dec using a 0?01 window (middle panels), and as a function 
of mean r PSF magnitude using a 0.1 mag window (lower panels). 
The mean difference is only calculated if there are at least 100 stars 
in the sliding window. The left hand panels correspond to proper 
motion in RA and the right hand panels correspond to proper mo- 
tion in Dec. We also plot the mean difference plus or minus the 
running 3cr-clipped standard deviation as the upper and lower sets 
of smaller points, respectively, in each panel. 

Figure|9]illustrates that the systematic differences between the 
HLC proper motions and those of GK04 are at a very small level, 
generally <2 mas yr~^. We note one clear systematic trend that 
the GK04 Dec proper motions are offset from the HLC Dec proper 
motions by < 2 mas yr^^, which is especially visible in the bottom 
right hand panel of Figure [9] We are currently unable to identify 
unambiguously the origin of the small systematic offset. The scat- 
ter in the proper motion differences (represented by the upper and 
lower sets of smaller points in each panel) is consistent with the 
stated proper motion uncertainties of ~3.9 mas yr^^ in GK04, and 
< 5 mas yr~^ for these particular stars in the HLC. 

Each light-motion curve in the LMCC has a different tem- 
poral coverage and number of epochs, a situation which is high- 
lighted in Figure |2] Consequently the uncertainties on the HLC 
proper motions exhibit a very inhomogeneous spatial distribution, 
and selecting proper motion objects using cuts on proper motion 
uncertainty results in a very inhomogeneous sample of objects. 
In Figure [To] we present the mean proper motion uncertainty in 
the HLC as a function of mean r PSF magnitude, number of 
epochs and time span of a light-motion curve for PSF-like objects 
(MEAN.OB JECT.TYPE = 6). The three panels from left to right 
correspond to different light-motion curve time spans of 3, 5 and 7 



years, respectively, while in each panel three curves are plotted, 
dotted, dashed and continuous, corresponding to <10, ~20 and 
>30 epochs, respectively. The curves in each panel show the mean 
proper motion uncertainty in the HLC as a function of mean r PSF 
magnitude. It is worth noting that ~76% of objects in the LMCC 
have at least 20 epochs, and that ~37% have at least 20 epochs with 
a time span of greater than 4 years. We find that 312819 objects in 
the LMCC (or ~8%) have proper motions /i with t-i/a-^ > 5. 



5 CONCLUSIONS 

The Light-Motion Curve Catalogue (LMCC) contains almost 
4 million light-motion curves for stars and galaxies covering 
~249 deg'^ in the SDSS Stripe 82. A light-motion curve provides 
wave band and time dependent photometric and astrometric quan- 
tities, where ~76 per cent of light-motion curves in the LMCC 
have at least 20 epochs of measurements. The LMCC is complete 
to magnitude 21.5 in r and i, and to magnitude 20.5 in mak- 
ing it the deepest large-area catalogue of its kind. The photometric 
RMS accuracy for stars is ~20 mmag at r ~ 1 8 mag and for galax- 
ies it is ~30 mmag at r ~18 mag. In both the RA and Dec coordi- 
nates, an RMS accuracy of ~32 mas and ~35 mas at r ~18 mag 
is achieved for pre-2005 data for stars and galaxies, respectively, 
and an RMS accuracy of ~35 mas and ^--^46 mas at r ~18 mag is 
achieved for 2005 data for stars and galaxies, respectively. 

The LMCC is thus an ideal tool for studying the variable sky, 
and in order to aid in its exploitation, we have created the Higher- 
Level Catalogue (HLC). The HLC consists of 229 derived photo- 
metric and astrometric quantities for each light-motion curve, and it 
is stored in a set of FITS binary tables, a format that is widely used 
by the astronomical community. The photometry presented in the 
HLC is fully consistent with the IV07 standard star catalogue and, 
since it is based on many more epochs, the random uncertainties in 
the mean magnitudes are smaller. Reassuringly, we show that only 
a very small percentage of standard stars (~0.065% or ~650 ob- 
jects) from IV07 are actually photometrically variable. The HLC 
proper motions of 30546 stars are consistent to wi thin uncertain- 
ties with those derived bv lGould & KoUmeieJ j2004h by combining 
SDSS DRl and USNO-B proper motions. 

The power in us ing the HLC is well illustrated by the work 
of IVidrih et all ilOOH) who construct a reduced proper motion dia- 
gram for Stripe 82 in order to find rare white dwarf populations, 
and consequently they identify 8 new candidate ultraco ol white 
dwarfs and 10 new candidate halo white dwarfs. Also, Becker et al.l 
|2008) report the discovery of an eclipsing M-dwarf binary system 
2MASS J01542930+0053266, having employed the corresponding 
light-motion curve along with radial velocity measurements to de- 
termine the masses and radii of the stellar components. We there- 
fore encourage the astronomical community to actively use and ex- 
plore these catalogues that are so ample in their content. 
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Figure 10. Plot of mean proper motion uncertainty in tlie HLC as a function of mean r PSF magnitude for PSF-like objects (stars). Tlie left, middle and right 
hand panels coirespond to light-motion curve time spans of 3, 5, and 7 years, respectively. In each panel, the dotted, dashed and continuous curves con'espond 
to <10, ss20 and >30 epochs, respectively. 
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